Level crossings in a cavity QED model 
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In this paper I study the dynamics of a two-level atom interacting with a standing wave field. 
When the atom is subjected to a weak linear force, the problem can be turned into a time dependent 
one, and the evolution is understood from the band structure of the spectrum. The presence of level 
crossings in the spectrum gives rise to Bloch oscillations of the atomic motion. Here I investigate 
the effects of the atom-field detuning parameter. A variety of different level crossings are obtained 
by changing the magnitude of the detuning, and the behaviour of the atomic motion is strongly 
affected due to this. I also consider the situation in which the detuning is oscillating in time and its 
impact on the atomic motion. Wave packet simulations of the full problem are treated numerically 
and the results are compared with analytical solutions given by the standard Landau-Zener and the 
three-level Landau-Zener models. 

PACS numbers: 45.50.Ct, 42.50.Pq, 42.50.Vk 
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I. INTRODUCTION 



The general theory of systems evolving with a pe- 
riodic Hamiltonians has long been known and under- 
stood. However, it is only recently that some of the 
theoretical predictions have been tested experimen- 
tally. It has turned out that the physical systems 
one first had in mind, namely these of solid state 
materials, have several drawbacks making it hard to 
realize such experiments. Instead of using the pe- 
riodic structure in crystal lattices, one may use a 
standing wave light field as an 'effective' potential, 
acting upon the internal electron states of an atom. 
Here the wavelength of the field sets the period, and 
due to lack of lattice defects in the standing wave and 
long atomic decay times, the decoherence times ex- 
ceed those in solid state systems by several orders of 
magnitude. One purely quantum mechanical effect 
of the atomic motion is seen when the atom is sub- 
jected to a constant force F. The optical potential 
then has a washboard shape. As long as the force 
is weak compared to the amplitude of the stand- 
ing wave, the quantum numbers that describe the 
system are still good numbers, apart from the fact 
that the quasi momentum will follow the classical 
evolution k = ko — Ft. Thus the system remains 
in its initial energy band while it sweeps the quasi 
momentum. Since the energy bands are periodic in 
k, the motion will also be periodic contrary to the 
constantly accelerating one expected from classical 
mechanics. This phenomena is called Bloch oscil- 
lations after the investigations 0, These oscil- 
lations have been verified experimentally in various 
systems; atoms in optical potentials |3- C3- 01 ■ BEC's 
in optical lattices |fj, |jj and solid state superlatti- 
cies yl. Some theoretical studies are presented in 
refs. |3, [^} , which analyze the problem from a Flo- 



quet or Bloch theory point of view. The dynamics 
can also be understood using a Wannier-Stark lad- 
der spectrum, consis ting of equally spaced discrete 
complex eigenvalues 

With no atom-field coupling, the atom moves 
freely and the spectrum, consisting of an infinite 
set of bare energy dispersion curves, is continuous 
E > 0. When the periodic potential is turned on, 
the degeneracy points where the curves cross are 
split and gaps are formed. The dispersion curves 
become periodic in k, with a period equal the Bril- 
louin zone. The avoided crossings are central for the 
evolution when k grows linearly in t. By lineare- 
lizing the dispersion curves around a crossing, the 
system serves as a testing ground for the Landau- 
Zener model ^3 El- This model treats a linear 
time-dependent two-level crossing system analyti- 
cally, and the asymptotic solutions for transitions 
between the two states are given. The Landau-Zener 
tunneling has been tested experimentally in connec- 
tion with Bloch oscillations 0, 0, 0] and also been 
the sub ject of theoretical investigations in the same 
system [l0| . 

In all the works mentioned above, the particle 
lacks internal structure. For example, in the situa- 
tion of atoms interacting with a standing wave field, 
it is assumed that the detuning between the atomic 
transition frequency involved and the field frequency 
is large, implying that one of the two atomic levels 
can be omitted after an adiabatic elimination. This 
results in a two photon process and no single pho- 
ton exchanges are present. The main reason for us- 
ing a large detuning is that only the ground state 
and not the upper atomic level is populated and the 
coherence times increase considerably. In the exper- 
iments, a strong driven light field is used making the 
decay of the field negligible. 

In this paper I investigate the situation when the 
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detuning may be small, which is fairly unanalyzed, 
one reference is 0| . It is understood that the paper 
is mostly of theoretical interest, neglecting decays of 
both the atom and the field. I do, however, men- 
tion how to overcome the problem of losses. The 
field may be driven by an external classical source, 
and one may use a A-type of atom where one of 
the lower states is coupled to the upper state by the 
standing wave and the other lower state is coupled to 
the upper one by an external classical traveling wave 
field. If the upper state can be adiabatically elim- 
inated, the effective two-level system has the same 
form as the one used in this paper, and thus the 
losses due to atomic excitations can be minimized. 
I argue how an effective model for describing losses 
of the field may be derived in certain case, and nu- 
merically show the effect of the losses. I keep the 
field quantized, assuming the field to be in an n 
photon Fock state; this does not make the analy- 
sis more complicated than for a classical field. In 
contrast to the paper I consider not only the 
zero detuning situation. By controlling the detun- 
ing, a new type of level crossing is achieved where 
three bare energy curves cross. This kind of crossing 
have similarities with the three- level Landau- Zener 
model 16], and comparisons between wave packet 
simulations of the full atom-field system and the an- 
alytically solvable three-level Landau-Zener model 
is carried out. Thus, this atom-field model provides 
one physically interesting system where this three- 
level Landau-Zener generalization is applicable. 

I also investigate the situation when the detun- 
ing becomes time-dependent and particularly when 
it oscillates in time, which is the situation when 
the atom is strongly Stark shifted by an alternat- 
ing external electric field. In these considerations, 
the force is not included, and hence the dependence 
of time arise from the oscillating detuning only. The 
atom again traverses a level crossing and is, in the 
adiabatic limit, transferred between different bare 
states, and, consequently, its motion may have an 
oscillating behaviour. 

The paper is organized as follows. In section ITT1 
the extended Jayness-Cumings Hamiltonian describ- 
ing the system without an external constant force is 
introduced. I discuss general features like the two 
natural basis states, the band structure and effec- 
tive parameters. The zero detuning case with the 
force applied is considered in lHIl and I show the re- 
sults from numerical wave packet simulations. The 
Landau-Zener transitions and the Bloch oscillations 
are discussed. The three-level crossing model is pre- 
sented in section IIVI and the wave packet simula- 
tions are compared with the analytic results of paper 
(l6| . In the next section I look at the situation 
when there is no force acting on the system, but 
the detuning is oscillating in time. It is found that 
a similar evolution may be obtained as in the case 



with a force. A discussion about how to minimize 
losses, of both the atom and the field subsystem is 
considered in section IVTI where I also show the ef- 
fect of field losses on the atomic inversion. Finally I 
conclude the paper with conclusions in lVIII 

II. THE MODEL 

The evolution of the combined atom plus a cavity 
mode is described by a Jaynes-Cummings type of in- 
teraction 0,0], where the coupling is x-dependent 
following from the variation of the electric field E(x) 
along the cavity mode, g{x) = d ■ E(x)/h. In partic- 
ular, the mode is assumed to have a standing wave 
shape, g(x) = 2g cos(qx). For cold atoms, their ki- 
netic energy term must be taken into account fully 
quantum mechanically |l9l |20| , and we arrive at 

h 2 d 2 A _ „ . , + , N 

2rhdi 2 ~ + J' Tz+9() 2cos (9 x H aV + aa )> 

(1) 

where m is the atomic mass, A is the atom-field 
mode detuning, (a) are raising (lowering) boson 
operators for the mode state, q = 2tt/X is the pho- 
ton wave number and <r z , are the Pauli z, rais- 
ing/lowering operators acting on the two-level atom. 
The dynamics governed by CI has been studied in 
numerous papers, see i.e [HlH I2I I2I I2I I2I l27j . 

In the following we will use scaled units. The pho- 
ton recoil energy = h 2 q 2 /2m defines a charac- 
teristic energy scale and q _1 sets a length scale. The 
new scaled parameters are 

x = qx p = 2 t = f^f 

(2) 



where the tilde ~ indicates unsealed quantities. 

The scaled Hamiltonian has a spatial period d = 
2ir, and according to Bloch or Floquet theory, it is 
known that the system eigenstates are determined 
by two quantum numbers; a discrete number v and 
a continuous number k [23 ■ The first v is the band 
index ranging from 1,2,... and k is the quasi mo- 
mentum taking values within the first Brillouin zone 
— 1 < k < 1 (scaled units). In some papers, the Bril- 
louin zone has been defined as —1/2 < k < 1/2, 
with the consequence that two sets of energ y b ands 
are obtained, one for each internal state |±) [lal2l|. 
The spectrum E u {k) forms an infinite set of continu- 
ous bands separated by gaps. For smooth potentials, 
the gap size decreases for increasing values of v\ in 
general the band-to-gap ratio falls off faster than 
\jv. Since upper | \) and lower | |) electronic states 
of the atom couple through absorption or emission 
of one photon, the number of excitations of the sys- 
tem is preserved by the Jaynes-Cummings type of 



3 



interaction JIJ. We may thus, for a given number of 
excitation, define the excitation eigenstates 



where \n) is the n photon number state of the mode. 
These will also be referred to as internal states. In 
each exchange of energy between the field and the 
atom, the atomic momentum is shifted by ±1 (in 
scaled units) from the mechanical 'kick' due to ab- 
sorption or emission of a photon. Consequently, an 
internal state, |+) (|— )), will only return to the same 
internal state after having shifted the momentum by 
an even number of momentum kicks. A state with 
momentum p — k, \k), and internal state |— ), is 
coupled to the following set of states 

~ \ \k + fi)\+), p odd, W 

where [i runs over all integers. Note that the mo- 
mentum of \ip^(k)) is k + fj,, and the quasi momen- 
tum k will therefore be imagined within the first 
Brillouin zone. There is an orthogonal set of states 
with even/odd reversed, which is obtained by using 
the initial state These two sets do not cou- 

ple, and in the special case of zero detuning 5 = 
their spectra become identical. This degeneracy will, 
however, never be interesting for us, since we will al- 
ways assume our system to be in one of the two sets. 
The states of Eq. |@J are egenstates of the Hamilto- 
nian in the absence of coupling go = 0, and I will 
refer to them as bare states. For a non-zero cou- 
pling the bare states are no longer eigenstates of the 
Hamiltonian, and the new eigenstates, 

H\<j> v {k)) = E v {k)\<f> v {k)), i/= 1,2,3,..., (5) 

will be called dressed states. In the theory of periodic 
operators, the dressed states are also known as Bloch 
or Floquet states and their properties are familiar. 
It is often convenient to transform between bare and 
dressed states 

l<M*0> = E M cHWM*)) 

(6) 

\MQ) = E„ c»(k)\Mk)), 

where the coefficients c^(fc) depend on k. 

The Hamiltonian has close similarities to the 
Mathiue equation [28|, and in the 5 — limit our 
'two-level' Hamiltonian is, actually, separable by a 
unitary transformation into two uncoupled Math- 
iue equations. In the opposite limit, go/ 5 — ► 0, 
the Hamiltonian is also diagonal in its internal 'two- 
level' structure. For the two disconnected equations 
obtained, the 'potentials' go as ~ cos 2 (a;), see 0. 
The 'potentials' equal | [1 + cos(2a;)], where the first 
term describe the situation when the two photon 
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FIG. 1: The lowest energy bands of the Hamiltonian 
Q as function of quasi momentum k for two different 
detunings S. The quasi momentum k runs between - 
1 and 1 in the plots. The Bragg scattering crossings 
are marked with circles and the Doppleron with square 
boxes. The coupling is Vo = 0.2 in both plots. 



kicks have opposite momenta, while the second term 
gives the processes when the absorbed and emitted 
photons kick the atom in the same direction, either 
left or right. In Fig. 1 two examples of the lowest ly- 
ing bands of the spectrum E v (k) for — 1 < k < 1 
are shown. In the first plot (a) the detuning is 
zero, 5 = 0, while the second plot (b) has a pos- 
itive detuning 5 = 0.7. The diagonal matrix ele- 
ments of the Hamiltonian, in the bare state basis, 
are s t _ l (k) = (k + ± 5/2, which are, of course, 
eigenenergies for zero coupling, and therefore called 
bare eigenenergies. For a non-zero coupling, the de- 
generacies are split, forming the gaps in the spec- 
trum. These avoided crossings will be of importance 
in the remaining of the paper. In the two plots we 
see two different kinds of crossings (there exist also a 
third type discussed in section IV) , those when both 
of the bare crossing energies have the same internal 
state |±) and those when they have different inter- 
nal states. In the first case we say that the crossing 
point is a Bragg resonance and the second is called a 
Doppleron resonance • Note that in (b) the most 
pronounced avoided crossing is between the second 
and the third band and not between the first and 
the second band. This is because the lowest bare 
eigenenergies that build up the two lowest energy 
bands both belong to states with the same internal 
state, and these states are only indirectly coupled by 
the Hamiltonian QJ. By controlling the detuning a 
varity of different spectra can be obtained. For ex- 
ample, the location of the gaps and also their size 
can be controlled to some degree. 

The spectrum or the dispersion curves E v (k), con- 
necting quasi momentum k with energy for a given 
band, contain information about the dynamics of 
the system. More precisely, the dispersion curves 
identify the behaviour of an atomic Gaussian wave 
packet in the presence of the periodic light field. For 
us, there are two types of wave packets that are of 
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interest; Gaussian dressed or bare state wave pack- 
ets. The first one is a Gaussian wave packet built 
out of dressed states according to 



dk tp v {k)\<t>v{k)) , 



where tp v {k) is Gaussian centered around some quasi 
momentum ko within the first Brillouin zone. Note 
that we have assumed only a single band to be oc- 
cupied, thus we do not sum over v. The second 
Gaussian state is 



I*) = / dp X (p)\p)\-) = 



E 



dkx(k + fi)\i^ f _ l (k)), 

and x(p) is a normalized Gaussian distribution with, 
for simplicity, a spread within the size of one Bril- 
louin zone. Using Eq. © for going between the two 
basis, it follows that \(k) = ^Pu(k — n)c£(k — fi). 
If c£J (k) is smoothly varying over the spread of the 
Gaussian this means that a Gaussian distribution in 
quasi momentum space gives a set of nearly Gaus- 
sian distributions, each shifted by unity in either di- 
rection, in real momentum space. The time evolu- 
tion of a dressed Gaussian wave packet is 

|$„(t)>= [ dk<p„(k)e- iE " ik)t \Mk)). (9) 



could, for example, be the force felt by a trapped 
atom due to gravity. Clearly, by including the ex- 
tra term, the periodicity is violated and the new 
eigenstates will have different properties than the 
former dressed states. The spectrum no longer con- 
sists of bands, but becomes purely continuous. It 
is also commonly known that the spectrum can be 
described by a discrete set of complex eigenvalues, 
called Wannier-Stark ladders . Contrary to what 
could be expected, the atomic wave packet will not 
constantly accelerate, but rather have an oscillatory 
behaviour. This characteristic has become known 
as Bloch oscillation originating from 0, , and it is 
usually understood either by Wannier-Stark ladders 
fll| or from the typical band structure of fig. 1 [jj. 
Here we will use the latter of the two approaches to 
analyze the dynamics. 

When the linear force F is weak enough, an ini- 
tial Gaussian dressed state will not populate nearby 
bands. It is the possible to show p| that the quasi 
momentum distribution obeys 



l\^k,t)\ 2 = -F^(k,t)\ 2 , (11) 



with solution 



\<p(k,t)\ a = \<p(k-Ft)\ s 



(12) 



By expanding the dispersion curve around ko and 
only including the first three terms (assuming E v (k) 
to vary little within the spread of tp v (k)), one may 
interpret the individual terms of the expansion: The 
first gives an overall phase-shift, the second and 
third ones are connected with group velocity and 
effective mass mi as 



dE v (k) 



dk 



k — kn 



and — = 

m 2 



d 2 E v (k) 



dk 2 



-ko 

( 10 ) 

In the absent of external forces, v g determines the 
velocity of the wave packet and mi determines how 
'fast' the wave packet spreads, see |53 for a detailed 
discussion. The wave packet thus behaves as freely 
evolving but with some effective parameters. This is, 
of course, all well known from the theory of electrons 
in crystal lattices. 



III. LANDAU-ZENER TRANSITIONS AND 
BLOCH OSCILLATIONS 

A very interesting phenomenon of the atomic mo- 
tion occurs when a linear force is allowed to act 
on the atom, i.e. the Hamiltonian (|TJ will include 
an extra term Fx, where the coefficient in scaled 
dimensionless units becomes F — FjqE^. This 



This is nothing but having an 'adiabatic' evolution. 
Within this limit, a Gaussian dressed state will move 
its mean according to k = kg — Ft, and as it exits one 
Brillouin zone it enters the same zone on the oppo- 
site side. All population remains in the same band 
as the initial state during its evolution, and we see 
from fig. 1 that, since the group velocity v g defined 
in eq. I|1(J|) will flip sign while traversing the quasi 
momentum, the atomic wave packet will oscillate in 
x-space. 

We know that the validity of adiabatic approxima- 
tions is related to the 'distance' between adiabatic 
energy eigen-curves |29| . Thus, the approximation 
that no other bands are populated will, most likely, 
break down near the avoided crossings. In this sec- 
tion we chose the zero detuning situation, S — 0. 
By neglecting coupling to other bands except close 
to a crossing, we may treat the problem as a two- 
level system with a lowest order effective Schrodinge 
equation 



8 



(k 



Ft 

Vo 



1 



2 



d„ 



where we have assumed 



(ko - Ft) 2 

(13) 

Doppleron resonance 



crossing. Here d M is the amplitude for the bare state 
IVvC^o — Ft)}- Thus, if we consider bands v and 
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v + 1, the various amplitudes are 

before crossing 



x 10 



^1/ 



C i/+1' 



after crossing 

before crossing 
after crossing 



(14) 



By linearalizing the energy dispersion curves around 
the first crossing and omitting the constant energy- 
shifts, we obtain 



d 







' Ft 


Vo " 
2 










v 
. i 


-Ft 




d^ 



(15) 



This is the analytically solvable Landau-Zener model 
[T3. fl3j . which has been widely used in different ar- 
eas of physics. The asymptotic solution at t = oo is 




600 1000 
t [h/E R ] 



(16) 



where the initial condition at t = — oo is |g^+i| = 1 
and A = ttVq /4F is the adiabaticity parameter. The 
breakdown of not having exact Landau-Zener transi- 
tions between neighboring states in the Bloch oscilla- 
tion scheme has been thoroughly studied in ^3] , and 
will not be the subject of this paper. The solution 
(|16fl is only an approximation, but, nevertheless, it 
gives a rough idea of the parameter dependence; a 
small force F, the velocity of the momentum wave 
packet is low, and a strong coupling Vo, the band gap 
is large and the energy curves are far apart, which 
gives an adiabatic transition between the two bare 
states |Vv(/«o + Ft)) and |^v + i(fc + Ft)). Thus, 
according to eq. <|14|) . no transition takes place be- 
tween the bands in such a situation. Population that 
escapes into higher bands due to non adiabatic evo- 
lution, will move more freely, since the higher bands 
are less coupled to one another. Note that in the 
Landau-Zener solution, the system is integrated be- 
tween — oo to +oo and it is therefore not fully accu- 
rate to identify the amplitudes \d^\ 2 with the ones 
in the Landau-Zener model. 

In fig. 2 we show the evolution of an initial Gaus- 
sian bare state in the presence of a force F. The 
atom is initially in its lower state with a spatial dis- 
tribution 



X(x) 



1 



(17) 



with the width = 50. Thus, it has an average 
initial quasi momentum fco = and in x-space it 
extends over several periods of the standing wave 
meaning that A| is much smaller than the Brillouin 
zone. Note that, since the initial state is a bare one, 
some of the higher bands will be populated already 
before the first crossing. However, the coupling is 



FIG. 2: Bloch oscillations of an initial Gaussian bare 
state in the lowest band for zero detuning 5 = 0. Con- 
tours with filled gray scales show the upper atomic wave 
packet and the non-filled contours the lower wave packet. 
The scaled coupling is rather weak in this example, 
Vo = 0.2, implying that bare and dressed states may have 
a large overlap. Thus, as the Gaussian quasi-momentum 
wave packet traverses the lowest band, the atom will 
Doppleron scatter at each Brillouin zone and therefore 
the internal state |±) is flipped, which is seen in the lower 
plot showing the atomic inversion (o" z ). Here the force 
F = 0.005. 



weak enough that only a few percentage of the pop- 
ulation is in higher dressed states. This is seen in the 
lower plot in fig. 2 where the atomic inversion (<j z ) 
is shown as function of the scaled time. From the 
atomic inversion we see that the atom 'jumps' be- 
tween upper and lower states and we also note that 
the evolution is very adiabatic and not much pop- 
ulation leaves the lowest band. As argued above, a 
larger force F should make the population of the 
oscillating part damp out. This is seen in fig. 3 
where we plot the same as in fig. 2 except that F 
is three times as large. Note that the Bloch period 
Tb = 1/F is a third as long, and how the atomic 
inversion damps out much faster than in fig. 2 due 
to the non-adiabatic transitions. 



IV. THREE-LEVEL LANDAU-ZENER 
TRANSITIONS AND BLOCH 
OSCILLATIONS 

The previous section dealt with the dynamics of 
a wave packet, evolving under the Hamiltonian Q 
with the additional 'force' potential Fx, in the sit- 
uation of zero detuning 5 = 0. Loosely speaking, as 
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400 
t [h/E R ] 



FIG. 3: This shows the evolution of the same initial 
state as in fig. 2 except that the force is three times 
as large, F — 0.015. The non-adiabatic transitions due 
to the increased force are seen both in the atomic wave 
packets and in the atomic inversion, which damps out at 
each transition. 



the wave packet approached a level crossing, it adia- 
batically absorbed or emitted a photon, resulting in 
a flip of internal state and a momentum 'kick'. We 
said that the atom was Doppleron scattered. The 
main difference between those Bloch oscillations and 
the ones most often studied is the internal structure 
of the system; with no internal structure, the adia- 
batic transition only results in two or no momentum 
'kicks' and no flip, which is that of Bragg scattering. 

The existence of Doppleron resonances in this 
model is just one example of how this model behaves 
differently from the standard systems lacking inter- 
nal structure of the particle. By varying the detun- 
ing, a third kind of resonance can be achieved, which 
is seen in fig. 4. In the figure 8 = 1 or 8 = — 1, and 
in both situations there is a level crossing between 
three bare energy curves. This scattering contains 
zero, one or two photon exchanges, giving a com- 
bination of both Bragg and Doppleron resonances. 
More generally we have a three level crossing for 
8 = (2j + 1), j = ... , —1, 0, 1, ... , where negative 
detunings are located at k — and positive ones at 
k = ±1. 

The evolution of an initial Gaussian bare state 
with the atom in its lower state when 8 = 1 is 
shown in fig. 5. From the atomic inversion, we again 
note that most of the population is within the low- 
est band. However, as the wave packet traverses 
through a crossing, the upper atomic state is be- 
ing populated. This does not mean that the sec- 
ond energy band is populated, but that close to a 




FIG. 4: The same energy bands as in fig. 1, but for two 
different detunings. For these detunings, 8 = 1, -1, there 
is a new kind of level crossing where three curves inter- 
sect. Hence, these crossings are neither purely Bragg nor 
Doppleron resonances. 



crossing the dressed state is a linear combination of 
bare states containing both upper and lower atomic 
states. The parameters are the same as in the ear- 
lier figs. 2 and 3, except that the force is weaker 
F = 0.0025. 



x 10 




1500 2000 2500 3000 
t [n/E tt ] 



FIG. 5: Bloch oscillations of an initial Gaussian bare 
state wave packet for 5 = 1 and ko = 0. Filled con- 
tour curves represent the atom in its lower state, while 
the non-filled curves (hardly seen in the plot) shows the 
upper atomic wave packet. From the lower plot, show- 
ing the atomic inversion, we note that the upper atomic 
state is populated almost only close to the curve cross- 
ings, see fig. 4. This does, however, not mean that the 
second band is populated when the quasi momentum 
wave packet passes a crossing, rather that the ground 
dressed state is a linear combination of both upper and 
lower atomic states close to a crossing. Other parameters 
are, F = 0.0025, V = 0.2 and A 2 X = 50. 

In order to get a deeper insight into the tran- 
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sition we may, similarly to the previous section, 
treat the Hamiltonian as an effective time depen- 
dent three level system by neglecting all other levels 
and putting k = fco — Ft. Let us shift the crossing 
to t = and the energy E — 0, and then we get the 
lincaralizcd Hamiltonian 



Hi = 



2Ft ^ 

v„ n vo 

2 " 2 

^ 

u 2 



2Ft 



(18) 



This is a three-level version of the Landau- Zener 
model and it has been solved analytically in a more 
generalized form in [lfij . I may introduce a tran- 
sition matrix, T, that takes the initial probability 
amplitudes squared at t = — oo to their final values 
at t = oo, provided that only one level is populated 
originally. In the above model, bare and dressed 
states become identical at t = ±oo, except for a flip 
in subscripts 1 <-> 3, and the matrix T has identical 
elements irrespective if it is given in bare or dressed 
basis. Thus we have for bare states 







- p2 


2P(1 - P) 
(1 - 2P) 2 


(1-P) 2 ~ 








2P(1 - P) 
_(1-P) 2 


2P(1 - P) 


X 


K + A 2 




2P(1 - P) 


P 2 





|d;-' 2 
Id 



oo |2 
ft+ll 



(19) 

where only one of the three levels is populated at 
time t — — oo. The probability is given by 



P 



(20) 



where the adiabaticity parameter is the same as in 
the original Landau-Zener model; A = 7rV r 2 /4P. 

A comment is in order about the result 1)19(1 . In 
the limit of t — > ±oo, the eigenvalues of the Hamilto- 
nian (|18|l are — oo,0,oo and the three states do not 
couple. However, with the full Hamiltonian Q, in 
the approximation k — fco — Ft, the crossings are 
periodic and the adiabatic eigenenergies are never 
infinitely far apart, so the states will always be cou- 
pled. If we start with an initial bare state wave 
packet in the middle between two crossings, corre- 
sponding to a time t = — t, the overlap with dressed 
states will be non-zero for several states. At t = — r, 
the state is propagated across the crossing till it 
reaches the next middle point at t = +r, and there 
the populations of different bare or dressed states 
are calculated. In the limit of weak coupling we may 
expect that only the three lowest bands will be pop- 
ulated after the passage of the crossing, assuming 
that the initial bare state corresponded to the low- 
est ground dressed state. These three probabilities 
are compared with the asymptotic ones \d'^L 1 „ M+1 1 2 
obtained from eq. I|19|) . One may expect that it 




FIG. 6: This figure shows the probability 12 U for the 
atom to be in its excited state |+) and with momentum 
— 1 < p < as function of force F (dotted line), after 
an initial Gaussian bare wave packet in the |+) state 
has been propagated across the crossing at k = seen 
in fig. 4. Dashed line displays the probability 1221 of 
ending up in the lowest band v = 1. The solid curve gives 
the result of eq. 1191 with initial condition \d~^l\ = 1. 
Finally the dotted-dashed line represents the result from 
numerical integration of eq. 11811 between t — — r and 
t = r corresponding to one half Brillouin zone (from 
fco = —1/2 to fco = 1/2). All four curves approaches 
zero when the force is increased, leading to a larger non- 
adiabatic transitions. The other parameters are Vo = 
0.2, fc = -0.5 and A? = 50. 



would be more correct to compare them with the cor- 
responding probabilities obtained from integrating 
the eq. (|18fl between — r and +r. The corrections 
of integrating over a finite time interval depends, 
of course, on how fast the solutions approache their 
asymptotic limits, for a thorough discussion see [Ioj|. 
In fig. 6, the dotted line gives the probability 

P+(-Kp<0) = J dp\(p,+\*(t = T))\ 2 , (21) 

where \^{t = r)) is the final state propagated across 
the crossing at k = shown in the right plot of 
fig. 4, as a function of the force F. The initial state 
is a Gaussian bare state with internal state |+) and 
fco = —1/2. Dashed line shows the probability of 
remaining in the lowest band 



P(v = l) = J ^dk\(<f> u = 1 (k)\V(t = T))\ 2 . (22) 

Finally the solid line and dotted-dashed curve dis- 
play IcJJJLjI 2 obtained with the transition matrix ac- 
cording to eq. (|19f) and from numerical integration 
of (I18|l between t = ±t respectively. It is seen that 
the probability from numerical integration follows 
the asymptotic solution very well. The coupling is 
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rather weak in this example, implying that bare and 
dressed states have a large overlap. 

I do not go into details in understanding the differ- 
ences between the three curves. In principle a similar 
analysis as the one performed in [lOj for the original 
Landau-Zener result could be done. However, the 
Dykhne-Davis-Pechukas method [3(], 0OJ for calcu- 
lating transition probabilities for two- level systems, 
which is used in |lCj is, obviously, not justified on 
our three-level system. 



V. CHIRPED ADIABATIC TRANSITIONS 

In the previous two sections the atom, while in- 
teracting with the standing wave field, was exposed 
to a constant force. The effect of this force could, 
in some linear or adiabatic regime, be effectively de- 
scribed by quasi momenta growing linearly in time. 
From the effective system Hamiltonian we could un- 
derstand how the dynamics of the atom behaves and 
how the transitions between states take place. In- 
stead of adding a force F to the system and con- 
sequently introducing a time dependence, we could 
externally vary either the coupling Vq or the detun- 
ing S during the interaction in order to insert explicit 
time-dependence. In this section we no longer apply 
an external force, but instead assume the detuning 
to be time dependent, and the change in time will be 
assumed slow, such that we remain in the adiabatic 
regime throughout the evolution. 




FIG. 7: The two lowest adiabatic energy bands E„(k, t), 
obtained by diagonalization of the Hamiltonian with a 
time dependent detuning S(t) = 5ocos(uit). It is seen 
that all crossings are avoided. The parameters are Vo = 
0.5, So = 2 and u = 0.1. 

One interesting choice of a time dependent detun- 
ing would be an oscillating one, S(t) — Socos(ujt), 
where w sets the time-scale. One could achieve 
this kind of detuning by, for example, Stark-shift 
the atom with an external alternating electric field. 



Usually the Stark-shift is quadratic in the field 
amplitude, so in order to have the above time- 
dependence the field variations must be chosen care- 
fully. By making an instantaneous diagonalization 
of the time-dependent Hamiltonian, we find time de- 
pendent adiabatic eigenvalues, E v (k,t). In fig. 7 we 
plot the time variation of the two lowest bands, for 
a fairly low detuning amplitude So = 2. For an adi- 
abatic process, a Gaussian dressed wave packet will 
stay within one band and evolve according to the ef- 
fective parameters (|1U|I . which will, of course, change 
in time. 

Assuming that we can approximate the system to 
contain just the two lowest dressed states (and k > 
0) we have the effective Hamiltonian 



(k- 



If- 

2 



2 2 



(23) 



and for simplicity we pick k 
Schrodinger equation 



d 



f2 









2 

Ml 

2 



1111 

2 



0.5, giving the 



(24) 





V2 




. Pi . 



The eigenvalues of the above Hamiltonian are 
±^^JS 2 (t) + Vq and thus, when S(t) changes sign 
we have an avoided level crossing. In the adiabatic 
limit, the population is transfered between the bare 
states across the crossing. The Schrodinger equation 
(|24|l with S(t) = Sq cos(uit) has been studied in [32| . 
The adiabatic solution, assuming an infinitely slow 
change, for the inversion is |3j| 



= 1 - 



1/2 



lV 2 + 5 2 (t)} 



1/2 



V$ cos (<Kt)) 



(25) 



si 



1/2 WS 



1/2 ' 



where 



dt' 



S 2 (t>) 



(26) 



and with 6(t) = Sq cos(wi), we can express this angle 
as an elliptical integral 



E(ojt\m) 
a/1 — m 



m 



m 



772-' (27) 



V 2 ' 
1 o 



where JH(u>t\m) is the elliptic integral of the second 
kind psj . In the adiabatic limit and first assuming 
cos (0(f)) = 1, the inversion H25|) will oscillate with 
period t osc = -k/uj. The oscillations coming from 
cos (4>(t)) will be seen as ripples on top of the other, 
slower oscillations, and the amplitudes of the two 
kinds of oscillations are determined by Vo and 5q. 
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As already mentioned, in the limit of large de- 
tuning, bare and dressed states are identical, and 
since the change is assumed to be adiabatic, the 
states will remain dressed also for small detunings. 
In fig. 8 we show the evolution of an initial bare 
Gaussian wave packet, when the detuning oscillates 
as 5(t) = Sq cos(wi). But since the initial detuning is 
rather large, the state has dressed features also when 
the detuning is small, which is seen in the oscilla- 
tions of the wave packet, see U^j. From the lower 
plot one notes that the oscillations of the inversion 
die out (compared to the stable adiabatic solution 
(J2SJ)), which is due to the non- adiabatic evolution of 
the wave packet. However, in terms of the transition 
rate between the two states, the adiabatic solution 
and the numerical one agree very well. The ampli- 
tude of the detuning is here S = 80, and by using 
typical experimental parameters (see next section) 
the unsealed amplitude is of the order of ~MHz. 
Achieving such strong Stark-shifts should be possi- 
ble 0. 




t [h/E R ] 



FIG. 8: The wave packet evolution for the time de- 
pendent detuning S(t) = So cos(u>i). It is seen that the 
wave packets moves back and forth as the detuning flips 
sign. This is due to the adiabatic transition between bare 
states with different group velocities v g . The oscillations 
in the wave packet are typical for Gaussian dressed states 
pi. H(i||. In the lower plot, the atomic inversion is shown 
from the numerical simulation (solid line) and for the 
adiabatic result i1'25ti (dotted line). The parameters are 
Vo = 10, So = 80, A x = 300 and u = 0.1. 



VI. THE PROBLEM OF LOSSES 

So far losses have been neglected, this is clearly a 
severe approximation, especially since the dynamics 
is assumed to be adiabatic. In this section I intend to 
discuss ways to minimize the losses, and also argue 
about the effects of losses of the field. 

Throughout the paper I have kept the field quan- 
tized, and assumed it to be in a Fock state. The 
Hamiltonian then has exactly the same form as it 
would have with a driven field with a large amplitude 
and the same mode shape. So the analysis is applica- 
ble to such a situation as well. One advantage with 
such fields is that losses are negligible. Thus, losses 
of the field may be decreased by drive the field with 
an external classical source. However, in order to 
have a stable periodic pattern of the field, it should 
be fixed using some mirrors/cavity. It is understood 
that since the upper atomic level is populated dur- 
ing the interaction it is most likely that losses will be 
of great importance. However we may estimate typ- 
ical time-scales in physical units for the system. In 
ref. p| the photon recoil energy is Eji — 1.03 x 10~ 10 
eV, which gives a characteristic time t — Tit/ « 1 
ms, where we took t = 200 from fig. 2. This rather 
long time may be decreased by using a stronger force 
F, which, however, increase the non-adiabatic con- 
tributions. As mentioned, losses of the field is mini- 
mized by using a strong driven field while one way to 
overcome the losses of the excited atomic state is to 
use a A-type of atom as in fig. 9, where the fields 
are assumed driven. When the detuning A in the 
figure is lar ge, the upper level can be adiabatically 
eliminated 36] to give the interaction Hamiltonian 
(scaled units) 



I Vq cos(a;) 
Vocos(x) — | 



(28) 



where the effective coupling is Vo — 0l ^ 223 . Thus, 
when the elimination of level |3) has been carried 
out, the dynamics of this model is exactly the same 
as for the two-level atom used in this paper. 

Finally one may use some arguments to get a 
rough idea about how the losses of the field affects 
the dynamical quantities. To solve the problem cor- 
rectly the master equation most be considered. How- 
ever, we may use some approximations to derive an 
effective model. By assuming that the field is ini- 
tially in a coherent state with a large amplitude, it 
is likely to assume that, in spite of the interaction 
with the atom, the field remains in some coherent 
state or in a linear combination of coherent states 
with the same amplitude. The effect of a zero tem- 
perature environment coupled to a coherent state is 
that the amplitude decrease as a — cto e ~ Kt ? where n 
is the coupling constant to the environment, see |37|. 
As the amplitude of the initial coherent field is large, 
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FIG. 9: A-configuration. The 2-3 transition couples 
to a standing wave mode, while the 1-3 transition to a 
traveling wave field. 



and hence the photon distribution is sharply peaked 
around its mean n, we introduce the effective cou- 



pling V n {t) = V£ 



with Vq 



VoVfi- Thus, the 



coupling is assumed to decrease in a similar way as 
the square root of the mean photon number n. This 
is, of course, a rough approximation; the evolution 
is still coherent, contrary to the one obtained from 
solving a master equation, and the spread of the 
photon distribution is not taken into account. How- 
ever, it should give a measure of typical decoherence 
times. In fig. 10 I plot the atomic inversion (a z ) 
as a function of time for three different couplings 
k _1 = 250, 500, 1000, where we have assumed an 
initial Fock state |1) of the field so no n dependence 
on the coupling. It is seen that, the larger coupling, 
the faster decay of the inversion; (cr z ) — > 0. This 
corresponds to typical decay times 1.5-7 ms. These 
are rather long for optical cavities (~ fis), but not 
for microwave cavities (~ 10 ms). In order to see the 
effect of the width of the photon distribution I have 
calculated (a z ) for various n (consequently various 
couplings, chosen as Voy/n/fie~ Kt ) and weight them 
with the distribution of a poisoning (coherent state). 
This has been done for the examples of fig. 10 with 
n = 50 and very similar results as in that figure was 
obtained, and are therefor not shown here. Hence, 
with initial coherent states with relatively large am- 
plitudes the same kind of evolution as for Fock states 
is expected. 



VII. CONCLUSION 

I have analyzed the dynamics of a single two-level 
atom inside a cavity, interacting with a standing 
wave mode field. By exposing the atom to a con- 
stant force of moderate amplitude it is possible to 
transform the time independent problem into one 
where the momentum grows linearly with time. This 
transformation is only valid for weak forces, sug- 
gesting that we can apply the adiabatic theorem, in 
other words, time can be considered as a parameter 



k =1000 




K =250 



200 400 600 800 

t [h/E R ] 



1000 



FIG. 10: This shows the results of having a decaying 
coupling Vq = V e~ Kt for n' 1 = 250, 500 and 1000. 
Here V e = 0.5, F = 0.02, <5 = and 5% = 50. The 
decay of the inversion is seen to increase for larger k, as 
expected. 



and the Hamiltonian can be instantaneously diag- 
onalized. The system, initially in an eigenstate of 
the Hamiltonian in the absence of the force, will re- 
main in the same energy band even when the force 
is affecting the dynamics; \<j) v {k = k — Ft)). As 
the quasi momentum is changed, the system tra- 
verses so called level crossings where the bare en- 
ergies cross. It is at these crossings that the adi- 
abaticity constrains are most likely to be violated, 
meaning that other states/bands may be populated. 
If the population still remains in a single band even 
through the passage of a level crossings, the atomic 
motion will be oscillating, contrary to standard be- 
haviour were the particle accelerates to higher and 
higher velocities due to the force. This is usually 
called Bloch oscillations and is a consequence of the 
energy gaps/crossings in the otherwice continuous 
spectrum. 

The Bloch oscillations of the atom are investigated 
for various values of the atom-field detuning. Bloch 
oscillations of atoms interacting with standing wave 
fields have long been known and also seen experi- 
mentally. But it has almost always been assumed 
that the detuning is large so that one internal elec- 
tronic state of the atom can be adiabatically elimi- 
nated, and the internal structure is overlooked. Here 
I have studied the situations when both atomic lev- 
els are included, and the differences in the dynamics 
due to the two-level structure. In the situation when 
just one level needs to be taken into account, the 
atom undergoes a two-photon process as it traverses 
a level crossing adiabatically, it is Bragg scattered. 
With small detunings, one-photon processes are also 
possible, making the atom flip its electronic state. In 
the zero detuning situation, we have seen that for a 
weak coupling and a small force the atom oscillates 
both in its motion but also between the two elec- 
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tronic states. We also note how non-adiabatic evo- 
lution started to dominate for larger values of the 
force. The Landau-Zener model was introduced in 
order to reach a deeper understanding of the param- 
eter dependences. 

For particular detunings a new kind of level cross- 
ing is present. In these situations, three bare en- 
ergy curves cross, and the atom may undergo zero, 
one or two photon exchanges with the field. The 
results from numerical wave packet propagations 
across such a level crossing has been compared with 
analytic solutions of the three-level Landau-Zcncr 
model. The agreement is very good even though 
there are numbers of approximations made in or- 
der to derive the effective three-level Landau-Zcncr 
Hamiltonian. 

In the section |3 we omitted the external force, 
and instead considered the atom-field detuning to 
be time-dependent 8 = 5(t). If the detuning changes 
sign, the bare energy curves, or in this adiabatic case 
surfaces, cross. By letting the initial Gaussian bare 
wave packet be centered around some momentum 
fco 7^ it has been possible to obtain a similar os- 
cillatory action as in the case with a constant force. 



The results are compared with the analytically ob- 
tained ones in the adiabatic limit. 

The effect of field losses and how to overcome such 
losses and also atomic losses were discussed in the 
last section IVII By using an A-type of atom and 
drive one transition by the cavity field and one by 
an external classical field an effective two-level atom 
with reduced losses can be derived by adiabatic elim- 
ination of the excited state. The field losses can be 
decreased by driving it with an external source. I 
display some numerical results when the field losses 
are included, and it is clear that the system is very 
sensitive to losses of realistic sizes. However, using 
the effective A model with driven fields, it does not 
seem too far out of reach to experimentally verify 
the theoretical predictions in this paper. 
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